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ABSTRACT 

We show that the extended main sequence turnoffs seen in intermediate age Large Magellanic Cloud 
(LMC) clusters, often attributed to age spreads of several hundred Myr, may be easily accounted for 
by variable stellar rotation in a coeval population. We compute synthetic photometry for grids of 
rotating stellar evolution models and interpolate them to produce isochrones at a variety of rotation 
rates and orientations. An extended main sequence turnoff naturally appears in color-magnitude 
diagrams at ages just under 1 Gyr, peaks in extent between and 1.5 Gyr, and gradually disappears 
by around 2 Gyr in age. We then fit our interpolated isochrones by eye to four LMC clusters with very 
extended main sequence turnoffs: NGC 1783, 1806, 1846, and 1987. In each case, stellar populations 
with a single age and metallicity can comfortably account for the observed extent of the turnoff region. 

The new stellar models predict almost no correlation of turnoff color with rotational vsini: the red 
edge of the turnoff is populated by a combination of slow rotators and edge-on rapid rotators. 

Subject headings: 


1. INTRODUCTION 

Exquisite photometry from the Hubble Space Tele¬ 
scope (HST) has recently revealed a surprise: many ^1- 
2 Gyr-old star clusters in the Large Magellanic Cloud 
(LMC) show extended or multiple mai n sequence turnoffs 


(MSTOs). F irst seen in NGC 1846 (Mackey & Broby 
Nielsen | 2007 ), and then also in NGC 1783 and 18(16 


] Mackey et al. 2008 Goudfrooij et al. 2009), these ex- 
tended turnoff regions now seem to b e a typical feature o f 
intermediate-aged LMC clusters (e.g. Milone et al. 2009). 
They cannot be explained away by photometric uncer- 
tainties or unresolved binaries, but can be reproduced 


by a spread of several hundred Myr in age (Goudfrooij 
etar][2MI| ). 


Large age dispersions are not expected in young and 
intermediate-aged globular clusters. They have generally 


not been seen in younger massive LMC clu sters (Bastian 
Niederhofer et al.|2015 ), thou^ Cor- 


& Silva-Villal2013 


renti et al. (|2'015 itound MCC 1856 to be consistent with 
an 80 Myr age spread at ^300 Myr. Young globular 
clusters have long been modeled with coeval stellar pop¬ 
ulation s of a single composition. |Elmegreen fc Efremov| 
(1997) showed that star formation in a typical ~10^ Mq 
cluster can finish in <10^ years, more than an order of 
magnitude smaller than the age spreads needed to ex¬ 
plain the MSTOs in the LMC. This picture of rapid star 
formation is confirmed by the h ierarchical distribution 


of clu ster ages seen in the LMC (Efremov & Elmegreen 
1998 ). Short bursts of star forination avoid the need 
to keep the cloud bound through the supernovae and 
copious ultraviolet radiation produced by the first gen¬ 
eration of stars. The existence of large age dispersions 
has also been challenged based o n the narrowness of the 
subgiant branch and red clump (Li et al. 2014 Bastian 
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fc Niederhofer||2015|), th ough this conclusion is disputed 
( |Coudtrooij et al.||2015 ). 

Stronger evidence tor multiple stellar populations has 
emerged for a few massive globular clusters. The clus¬ 
ter oj Cen has two well-separated track s in the color- 
magnitude diagram (Bedin et al. 2004), which can be 
interpreted as a ve ry large helium en hancement in the 
smaller population ( Piotto et al.|2005 ). M2 also has sev¬ 
eral cl early separat^ tra cks in a color-magnitude dia¬ 
gram (Milone et al. 2015), each cha racterized by a dif - 
ferent chemical abundance pattern (Yong et al. 2014). 
Similar results have been fou nd for a handful of othe r 
Galactic globular clusters (e.g. Marino et al.|2008' |2009|). 


Many authors have explored rotation as an alterna- 
tive explanation for the extended MSTO in intermediate- 
aged clusters. Rotation has long been known to affect the 
MSTO due both to an inc r ease in main sequence lifetime 
(Meynet & Maeder 2000 Ekstrbm et al.j 2012) and to 
viewing angle differences ( von Zeipel|1924[ ). Extreme ro¬ 
tators appear to be ubiquitous among early-type stars 
in the Solar neighborhood. Ve ga (a Lyrae) is observed 


to rotate at ^90% o f breakup (Aufdenberg et al. 2006_ 
Peterson et aH2006 ); Altair, a Oph (Kasalhaguej, and a 


Cep (Alderamin), three of the brightest, closest A-type 
stars in the sky , are all similarly rapid rotators (Mon- 
nier et al. 20071 Zhao et al. [2009). These stars are sig- 


niticantly deformed, with oblatenesses of ^0.8-0.85 and 
large pole-equator differences in temperature. 

The effect of orientation on the M STO in intermediate- 
age LM C clusters was explored by [Bastian ik de Mink] 
(2009|), who found that it could explain at least some of 
the observed spread. The ability of rotati on to explain 
the entire spread was directly challenged bylGirardi et al.| 

([MTl). 


Girardi et al. (2011) and 
both attempted to model LMC c 


Yang et al. (2013) have 

usters with coeval pop¬ 


ulations of rotating stars. Girardi et al. (2011) found 
that two rotat ion rates could no t explain the observed 
turnoffs, while Yang et al. (2013) also found rotation to 
provide only a partial solution to the extended MSTO. 
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Goudfrooij et al. (2014) argued in favor of age rather 
than rotation on dynamical grounds, based on a corre¬ 
lation between turnoff width and central escape velocity 


a proxy for the c luster’s ability to retain gas). Goud¬ 
frooij et al. (2014) also found a correlation between the 
prevalence or the bluest turnoff stars and of secondary 
red clumps, as e xpected from an age disp ersion. 


More recently, Brandt & Huang (2015b) noted that the 
nearby Hyades and Fraesepe opeii clusters seem to be in¬ 
compatible with a single nonrotating stellar population. 
Stars at the top of the main sequence turnoff appear, 
at more than 99.99% significance, to be younger than 


new rotating stellar models of 

Georg 

V et al. (2013b) as 

described in 

Brandt & Huang 

2015a 

), Brandt &: Huang 


(2015b) were able to ht both clusters with a single stel- 


lar population at a range of initial rotation rates and 
orientations. We use similar techniques in this paper to 
investigate the effect of rotation on the MSTO in much 
richer LMG clusters. 

We organize the paper as follows. In Section we 
review our method of interpolating the rotating stellar 
models and performing synthetic photometry. This pro¬ 
cedure is desc ribed in much more detail in IBrandt 
Huang (2015a), hereafter BH15. Sectioncompares our 
approach to previous attempts to explain the effect of 
rotation on the MSTO. We describe our results in Sec¬ 
tion 1^ and discuss their implications on the width of the 
subgiant branch and red clump, and on Galactic clusters, 
in Section We conclude with Section where we also 
suggest several areas for future research. 

2. METHODOLOGY 

Our methodology is described in detail in BH15. 
Briefly, we interpolate rotating stellar models and per¬ 
form synthetic photometry with a full gravity darkening 
calculation to obtain isochrones and color-magnitude di¬ 
agrams. We review the procedure here and refer the 

reader to BH15 for a more thorough discu ssion. _ 

We t ake the rotating stellar models of [Georgy et al. 


(2013b) as our basic inputs. We compute synthetic pho¬ 


tometry using the gravi ty darkening model of [Espinos^ 
Lara & Rieutord (2011). We integr ate the specific in¬ 
tensiti es of the model atmospheres ofjGastelli fc Kurucz 


( 12004), convolved wit h HST filter transmission curves 


(Siriahni et al. 2005), over the projected stellar sur- 
tace. We then interpolate the rotating models in mass 
and metallic ity using the fine nonrotating grid of jGirardij 
et al. (2002). We compute a rotational correction factor 
for each color and for each rotating model, and interpo¬ 
late these factors linearly in magnitude (as power laws 
in flux) between the rotating models. We then apply 
the interpolated rotational correction factors to nonro¬ 
tating models at the same points in their evolution and 
scale the main sequence lifetimes to match those of the 
rotating models. This procedure preserves the exact ro¬ 
tating models at the masses, metallicity, and rotation 
rates where they exist, and simply interpolates between 
them. 


The Georgy et al. (2013b) rotating models are only 
available from 1.7 to 15 Mq; they resolve neither the 
transition from a convective to a radiative core, nor the 
transition from a radiative to convective envelope. Both 
transitions happen below 1.7 Mq and have a large impact 


on stellar rotation. A convective core is uniform in com¬ 
position, so that mixing at the outer core boundary effi¬ 
ciently supplies additional nuclear fuel to the central re¬ 
gion where it is burned. Stars with convective envelopes 
shed angular momentum in a magnetized wind and spin 
down rapidly. Without the models to resolve this region 
of parameter space, we simply extrapolate the rotational 
correction to nonrotating models down to 1.45 Mq. Our 
results at ^2 Gyr ages, where the turnoff stars are ~1.5 
Mq, should therefore be viewed with caution. 

BH15 constrained the age and composition of coeval 
star clusters by computing and multiplying the poste¬ 
rior probability distributions of individual stars. This 
method neglects binaries, which could be observation- 
ally excluded from nearby clusters. The clusters we study 
here are much too distant to resolve binaries either inter- 
ferometrically or spectroscopically. We do not attempt 
to generalize the method of BH15 to statistically account 
for unresolved binaries, deferring that to a future paper. 

3. GOMPARISON TO PREVIOUS STUDIES 


Several papers have previously attempted to explain 
the extended MSTOs seen in LMC cluste rs by a distri¬ 
bution of rotation rates and orientations. IBastian fc del 


Mink (2009) suggested that the rapid rotation of massive 
stars could modify the morpholo gy of isochrones and lead 
to an extended MSTO. However, Bastian fc de Mirikj only 
considered orientation effects and small differences in ef¬ 
fective temperature, ignoring the substantial increase in 
main sequence lifetime that rotation can provide, and 


their results were directly challenged by Girardi et al. 

(IMIt. 


' Girardi et al. (2011) performed stellar evolution calcu¬ 
lations with two rotation rates, 0 and 150 kms“^, and 
without core convective overshoot. They also computed 
nonrotating models with a core overshooting parameter 
of 0.25. They found that the combination of rotating 
and nonrotating tracks could not reproduce the extended 
MSTO of the LMG cluster NGG 1846. Only a very wide 
distribution of ages combined with core ov ershooting pro¬ 
vided a satisfactory fit. The results of [Girardi et al.j 
([2011[) contradict some of the approximations made by 
Bastian h de Mink (2009), who did not run full stellar 

evolution calcula tions. _ 

More recently, Yang et al. (2013) investigated the ef¬ 
fects of rotation on the iVlSl'U using their own stellar 
evolution tra cks computed from the Y ale Rotation Evo¬ 
lution Code ([Pinsonneault et al.|1989[). They found that 
stellar rotation could lead to extended MSTOs in clus¬ 
ters with ages between 0.8 and 2.2 Gyr, but could not 
account for the full width of the MSTO in NGG 1806. O 


et al. (2015) and Li et al. (2015) attempted to reproduce 
the color-iriagnitude diagram of several clusters with a 
range of ages, binary fractions, and rotation ra tes; their 
best-fi t models retained large spreads in age. [Li et al. 


genera ted two evolutionary tracks from Georgy 


et al.j (2013b) to explore the effect of rotation on post- 


mam-sequence stars. They concluded that the effect of 
rotation on the morphology and luminosities of the sub¬ 
giant branch should be negligible due to the expansion 
of the stellar envelope. 

In this work, we use the stellar evolution models of 


Georgy et al.j (2013b), taking the extended main sequence 
litetimes and orientation effects fully into account. Un- 
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like the mode ls used by 
et al. (2013b) tracks ina 


Girardi et al. (2011), the Georgy 
ude modest core convective over- 
shoot (a = 0.1) in both the rotating and nonrotating 
cases, and cover a range of initial rotation rates. We 
make no attempt to model the post-main-sequence por¬ 
tion of the color-magnitude diagram. _ 

Our approach is similar to that taken by |Yang et al. | 
(2013), though we adopt dif ferent stellar evolu tio n cal¬ 


culations (se e Section 2.3 of Yang et al.| (2013)). Yang 
et al. (2013) also only considered initial rotation peri- 


ods of 0.49 days except in the case of NGG 1806, for 
which they assumed a Gaussian distribution centered 


on On /Orri t. = 0.3 with a dispersion of 0. 06. Li et al. 


( 2015 ) and Li et al. (2015), like us, used the Georgy et al. 


(2U13b) models. However, they only applied a rotational 


correction to other stellar evolution models with a fixed 
distribution of rotation rates, and do not appear to have 
accounted for gravity darkening and orientation effects. 


4. RESULTS 

4.1. The Main Sequence Turnoff With Rotation 

While ro tation has long been know n to prolong stellar 
evolution (Meynet & Maeder 2000), large sets of stel- 
lar evolution models have only recently bec ome available 
(Ekstrdm et al. 2012 Georgy et al.|2013a|bl). These mod- 


els allow us to quantity the spread ot the MSTO due to 
rotationaUy enhanced main sequ ence lifetimes. 

Eigure fl] shows the results of Ekstrom et al. (2012), 
who computed tracks with initial rotation rates 0 and 
^58% of critical, in physical units of Teff and log L (left 
panel) and in synthetic HST photometry at cost = 0.1 
(i « 26°, right panel). The main sequence turnoff has 
a pronounced extent in all populations younger than ~2 
Gyr. At older ages, the turnoff stars had outer convective 
zones while on the main sequence and shed their angular 
momentum while young. 

The left panel of Eigure is not observable. Golors 
are nonlinear functions of 5fective temperature and of 
metallicity; a turnoff region that is extended in physical 
units can be much less extended in photometry. There 
is no reason for the extent of the main MSTO to remain 
roughly constant below ^1.6 Gyr as it does in the left 
panel of Figure Indeed, the right panel of Figure 
shows a much less extended MSTO at ages <500 Myr. 

Figure shows the area of the MSTO on a color- 
magnitude diagram, with two choices of HST filter sets 
and with units of mag^ (magnitude times color), at 
and Z = 0.008. We define the turnoff area here as 
the area enclosed by a nonrotating isochrone and an 
isochrone rotating at Uo/Ocrit = 0.6 (Oo/flcrit = 0.58 
for Z = Zq = 0.014) viewed at cosi = 0.9 [i = 26°). 
An initial rotation rate 60% of critical appears to be rel- 
ativel y common for ~2 stars in the Solar neighbor¬ 
hood ( Zorec fc RoyerpOlY ). As noted in Section[^ Vega, 
Altair, a Gep, and a Uph all rotate much faster than 
this. We show two HST filter combinations: F435IU 
and F8UW (solid blue line) and F555IU and F814IU 
(dot-dashed red line). The three filters FAZbW, F555IU, 
and F814IU correspond roughly to Johnson R, U, and 
/, respectively. 

The MSTO in the physical Hertzsprung-Russell dia¬ 
gram (left panel of Figure]^ is roughly constant in area 
on the log-log plot for ages <1.6 Gyr. In the observable 


units of Figure the extent of the turnoff region has 
a clear peak between ^I and 1.5 Gyr, and is much less 
extended at both younger and older ages. There is no 
extended turnoff at young ages because B — I and V — I 
depend weakly on Tefj at high temperatures, and espe¬ 
cially at low metallicities. The metallicity dependence 
means that the extended MSTO phenomenon should ap¬ 
pear at a narrower range of ages in the LMG than in 
the more metal-rich Milky Way. At old ages, the turnoff 
stars shed their initial angular momentum early in their 
main sequence lives and become effectively nonrotating. 

We caution that the decline shown in the turnoff area 
around 1.5-2 Gyr cannot be trusted in detail for the 
Z = 0.008 isochrones. The Zq models (open star sym¬ 
bols) clearly show this decline at ~2 Gyr, but the low 
metallicity models only extend to 1.7 Mq. We extrap¬ 
olate them to 1.45 Mq for our analysis, implicitly as¬ 
suming that the effects of rotation change little from 1.7 
to 1.45 Mq. Depending on the details of the transi¬ 
tion from a convective to a radiative core and from a 
radiative to convective envelope, the age at which the ex¬ 
tended MSTO disappears could shift, and its disappear¬ 
ance could become more or less abrupt. The Zq points, 
which are computed from full isochrones and are valid 
at all ages, show an even more dramatic falloff in MSTO 
area around 1.5-2 Gyr than the extrapolated Z = 0.008 
models. 


4.2. Application to LMC Clusters 

We now attempt to reproduce observed intermediate 
age LMG clusters with very extended main sequence 
turnoffs, the best evidence for extended star formation 
in young globular clusters, with coeval stellar popula- 
ti ons. We simp l y inter polate the rotating isochrones 
of Georgy et al. (2013b) as described in Section and 
perform our owh synthetic photometry. We adopt the 
rich clusters NGG 1783, NGG 1806, and NGG 1846 as 
our main comparison sample. These were three of the 
first intermediate-aged clusters to show multiple MSTOs 


(|Mackey & Broby Nielsen|2007| |Mackey et al.|2008|), and 
show some of the most dramatically extended turnoffs. 
We also show the you nger, bluer cluster NGG 1987, 
which was claimed by Milone et al. (2009) to have a 
higher metallicity, lower extinction, and larger fractional 
age spread than other LMG clusters. 

Figure overlays our interpolated isochrones on the 
decontaminated HST color-magnitude diagrams of NGG 
1783, 1806, 1846, a nd 1987 as reduced and cleaned by 


Milone et al. (2009). All isochrones have been shifted by 
hand using a reddening and distance modulus to bring 
them into rough agreement with the cluster observations. 
We have made no attempt at a detailed parameter search. 
The line colors represent initial rotation rates from zero 
to 90% of critical, while the line types show viewing an¬ 
gles from cosi = 0.1-0.9 (i « 26°-84°). The three view¬ 
ing angles are indistinguishable in the color-magnitude 
diagram for rotation rates Dq < 0.3Dcrit- 
We emphasize that an initial rotation rate Do/Dcrit ~ 
0.9 does not appear to be exceptional among early-type 
stars near the Sun. Vega, Altair, a Cep, and a Oph all 


0.01, 

0.5, 0.5, and 0.9, respectivelv (^ 

4ufdenberK et al. 

2006 

Monnier et al. 

2007 Zhao et al 

2009|. We have 
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Fig. 1.— The area bracketed by nonrotating and Ho/f^crit = 0-58 isochrones at Zq (|Ekstrom et aL|2012|l, in units of logL and Tgff (left 
panel) and in the HST filters E435VF and F814VF (roughly Johnson B and 7) at p = cos 2 = U.i (rigtit panel). The shaded area in physical 
units is relatively constant with age until it disappears around 2 Gyr. At this age, the turnoff stars have spent their main-sequence lives 
with convective outer envelopes, and those with rapid initial rotation have spun down. The area in photometric units is not constant with 
age due to the nonlinear relationship between color and Teff. The relation between color and Tgff also depends strongly on metallicity. 



Fig. 2.— Area bracketed by the nonrotating and r2o/f2crit = 0.6, 
cosi = 0.9 isochrones (Qo/f^crit = 0.58 for Z = Zq = 0.014), 
which we adopt as the extent of the MSTO. In both of the HST 
filter combinations shown, and particularly at low metallicity, the 
extent of the turnoff is very small below ~500 Myr and peaks 
between 1 and 1.5 Gyr. The details of the disappearance of the 
extended MSTO in the Z = 0.008 isochrones at ~1.5-2 Gyr should 
be vi ewed with caution. T he Zq isochrones (open star symbols, 
from |Ekstrom et al.||2012|) resolve the transition from convective 
cores and radiative envelopes to radiative cores and convective en¬ 
velopes, but the low-metallicity isochrones (solid and dot-dashed 
lines) do not. 


these stars’ values of C and i on Figure 

Table summarizes our adopted cluster parameters. 
All isochrones assume Z = 0.008, reddening of A435VF — 
FS\4W = 0.09 to 0.18, and true distance moduli of 
18.60 to 18.65 mag assuming Fitzpatrick (19991 dust with 
Rv = 3. This is larger by ^0.1-0.15 m ag than the best 
LMC distance (Pietrzyiiski et al.|2013 ), a discrepancy of 
^5-7% in linear distance. Such a disagreement could be 
from, e.g., systematic errors in our synthetic photometry, 
an incorrect assumed metallicity, or incorrect assumed 
dust properties. 

Isochrones with a variety of rotation rates and viewing 
angles, but with a single age and composition, can eas¬ 
ily reproduce the observed extent of the MSTO. These 
results stand in stark contrast to those of lOirardi et al.l 
(2011), who concluded that rotation could not explain 


TABLE 1 

Adopted Parameters for Four LMC Clusters 


NGC 

O 

1 

E(B-I)* 

Z 

Age (Gyr) 

1783 

18.65 

0.09 

0.008 

1.55 

1806 

18.65 

0.10 

0.008 

1.55 

1846 

18.65 

0.12 

0.008 

1.55 

1987 

18.60 

0.18 

0.008 

1.00 


* More precisely, i?(F435VP — ^81414^) 


the turnoff of NGC 1846. With a modest core over- 


shoot parameter of 0 .1 in all models (Ekstrom et al. 2012_ 
Georgy et al.|201^), rapidly rotatirig stars spend ~25% 


longer on the mam sequence than nonrotating stars. This 
accounts for most of the turnoff width and mimics a 
spread of ~25% in age (up to ~400 Myr at the age of 
these clusters). Orientation effects further increase the 
observed width of the turnoff region. We note that a 
fractional age spread of ~20-30% is roughly consistent 
with the claims for LMC clusters (80 Myr at a n age of 
300 Myr, yd several hundred Myr at 1-1.5 Gyr|Milone| 


et al. 2009 Rubele et al. 2011 2013 Correnti et al.|2015 ). 

Measurements of the projected rotational velocity 
V sin i co uld help to distinguish the effects of age and 


rotation. Bastian & de Mink (2009) identified the red 
end of the turnott with the most rapid rotators, pre¬ 
dicting a strong correlation of vsini with color. The 
models we present, however, show a much weaker corre¬ 
lation. Figure shows the initial angular momenta and 
current projected rotational velocities of one realization 
of a 1.55 Gyr-old, Z = 0.008 population assuming 10 
mmag (1%) photometric errors. We use a Gaussian dis¬ 
tribution of rio/flcrit centered on 0.5 with a dispersion 
of 0.3, in ro ugh agreement with t he values for nearby A 
and B stars (Zorec & Royer 2012). We neg lect the effect 


of gra vity darkening on the Interred v sin i (Fremat et al. 
200 ^. 


Unlike Bastian & de Mink ([2009 ), we find the red end 
of the turnott to be populated by a combination of slow 
rotators and edge-on rapid rotators, while the blue edge 
is comprised of pole-on rapid rotators. The dispersion of 
V sin i increases with redder colors, but the correlation of 
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Fig. 3. — Color-magnitude diagrams of four intermediate-age LMC clusters with strikingly extended MSTOs: NGC 1783 (top left), 1806 
(top right), 1846 (bottom left), and NGC 1987 (bottom righC. The various curves are all coeval Z = 0.008 isochrones at different initial 
rotation rates (colors) and viewing angles (line types). Table[l]lists the assumed cluster ages, reddenings, and distances. We have indicated 
bright, nearby A stars with Q and i corresponding to the most rapidly rotating, flo/flcrit = 0.9 isochrones; all three are among the 60 
brightest stars in the sky, and the 12 brightest A stars. A range of stellar rotation rates and viewing angles can comfortably account for 
an extended MSTO at a single age and metallicity. 


color with mean v sin i is weak. Intriguingly, this agrees 
with the weak correlation of v sin i with color seen in the 


.^1.6 Gyr-old Galactic cluster Trumpler 20 (Platais et al. 


20121. The difference in mean B — I color between stars 
with tisini < 100 kms“^ and those with tisint > 100 
kms“^ in Figure]^ is just 2 mmag. Additional realiza¬ 
tions of the same cluster show that these B — I mean 
color differences are consistently <10 mmag. As the left 
panel of Figure]^ shows, even the initial angular momen¬ 
tum Go/Gcrit is only moderately correlated with MSTO 
position. 


5. DISCUSSION 


A major objection to age spreads in LMC clusters is 
the narrow distr ibution of stars along the post-main- 
seque nce tracks (Li et al. 2014 Bastian fc Niederhoferl 


2015 1, though this result is disputed 


20151. Rotating stellar models extend the IVIS'LO region 


et al. 


in two ways: by increasing the main sequence lifetime by 
a rotation-dependent factor, and by adding variation in 
color and apparent luminosity with viewing angle. The 
former effect introduces a range of turnoff masses at fixed 
age, and thus a range of stellar masses on the subgiant 
branch and in the red clump. Such a range of masses is 
exactly what would result from a spread in ages rather 
than rotation rates. 
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Fig. 4.— Initial angular momenta f^o/^crit (left panel) and current projected rotational velocities vsini (right panel) for one realization 
of an unreddened, 1.55 Gyr-old, Z = 0.008 stellar population with 10 mmag photometric errors in all bands. We assume a Salpeter IMF and 
a Gaussian distribution of flo/^crit with a mean of 0.5 and a dispersion of 0.3 (truncated so that flo/f^crit ^ (0? 1)) and neglect binaries. 
The re d side of the turnoff is populated by a combination of slow rotators and edge-on rapid rotators. In contrast to |Bastian &: de Mink| 
l |2009| |, who identified red colors with rapid rotation, the models we present show no significant trend of vsini with color, 'i'he ditierence 
in mean B — I color between stars with vsini < 100 kms“^ and those with usini > 100 kms“^ is consistently <10 mmag for realizations 
of this cluster. 


Figure shows that much of the extended MSTO phe¬ 
nomenon in our scenario arises because of the increase 
in lifetime with rotation (the line colors indicate dif¬ 
ferent initial rotation rates). Accounting for a broad 
MSTO with a spread in rotation rather than age can 
only reduce the spread in turnoff masses by a modest fac¬ 
tor, corresponding roughly to the fraction of the spread 
attributable to orientation. The difference in turnoff 
mass between a nonrotating and a very rapidly rotat¬ 
ing isochrone is nearly 10%, similar to the mass range 
produced by a 25% age spread. 

Another puzzle is the apparent lack of extended 
MSTOs in rich, ~1.5 Gyr-old Galactic open cl usters. 
NGC 7789 has a slightly super-Solar metallicity (Over- 
beek et al. 2 015) and an estimated age of 1.4-1.6 Gyr 
(|Kalirai et al.||2008 Gim et al.||1998 1, but does not have 
an obviously extended MS'l'U. 'I'rumpler 20, like NGG 
7789, has a slightly super-Solar me tallicity dGarraro et al . 
2010) and an age ^1.4-1.7 Gyr (Carraro et al.||2014p. 
Trurnpler 20 does show hints of an extend ed MS'l (J cor¬ 
respo nding to a ^300 Myr spread in age (Garraro et al. 
2014), tho ugh much of this m ay be due to ditterential 


reddening (Platais et al. 


20121 . 


The lack ot extended IVlS'l'Us in Galactic clusters could 
simply be due to the fact that the turnoff width is sharply 
declining at ages of 1.5-1.8 Gyr at Zq (Figure [^. NGG 
7789 and Trurnpler 20 have super-Solar metmlicities, 
and could plausibly show a decline in turnoff width at 
younger ages (equivalently, at lower turnoff masses). The 
rotating models we use also suggest that these clusters 
are slightly older than has been reported in the literature; 
the expected MSTOs would then be correspondingly nar¬ 
rower. 


One way to simultaneously account for a narrow sub¬ 
giant branch and for different turnoff behaviors in the 
Galaxy and the LMC is with systematically faster ro¬ 
tation at low metallicity. The opacity in stars and in 
molecular clouds is dominated by heavy elements. This 
opacity, in turn, sets the temperature gradient at which 
convection begins. Changing the heavy element abun¬ 
dance by a factor of three, roughly the difference be¬ 
tween young Galactic open clusters and the LMC, could 
plausibly have a strong effect on the dissipation of an¬ 
gular momentum. A bias toward rapid rotation at low 
Z would increase the relative importance of orientation 
in accounting for the turnoff width. It would thus allow 
for a narrower distribution of turnoff masses, while also 
narrowing the MSTO in higher-metallicity populations. 


6. CONCLUSIONS AND FUTURE WORK 

Our results show that variable stellar rotation and ori¬ 
entation at hxed age and composition can easily account 
for extended MSTOs, even in those LMC clusters with 
the most dramatically extended turnoffs. Stellar rotation 
naturally explains the appearance of a broad MSTO at 
ages just younge r than 1 Gyr and its disa ppearance at 
ages of ~2 Gyr. Brandt & Huang (2015b I showed that 
rotation can also remove a spread in color-magnitude 
ages in the much closer and less massive Hyades and 
Praesepe clusters. Taken together, these results strongly 
support the formation of open clusters in short bursts of 
star formation, forming essentially coeval stellar popula¬ 
tions of uniform composition. We suggest several areas 
of further study to validate and refine this picture. 

In principle, with a sufficiently fine grid of stellar mod¬ 
els and a careful accounting of stellar binaries, we can 


































7 


invert the observed distribution of stars in the color- 
magnitude diagram to obtain the distribution of initial 
rotation rates. We have made no attempt to do so in 
this paper. Such an inversion would require marginal¬ 
izing over parameters including reddening, distance, bi- 
narity, metallicity, age, and rotation, a highly complex 
process with strong covariances between the parameters. 
The ~0.1-0.15 mag difference between the best LMC dis¬ 
tance and the distance moduli used in Figure]^ though 
modest, suggests that we would need to treat any result 
with caution. 

Another promising avenue for future work is a detailed 
study of the disappearance of the extended MSTO phe¬ 
nomenon at ages ~2 Gyr. We are unable to res olve this 
transi t ion in t he low-metallicity case because the [Georgia 


et al. (2013b) rotating stellar models only extend to 1.7 
Mq. We extrapolate them down to 1.45 Mq, but in 
order to resolve the change in MSTO morphology, we 
would need models at many rotation rates extending at 
least to masses of ~ 1.3 Mq. Two important transitions, 
from a convective to radiative core and from a radiative 
to convective envelope, occur below 1.7 Mq. Both tran¬ 
sitions change the effect of rotation on stellar evolution. 
Our extrapolation, even down to 1.45 Mq, makes it dif¬ 
ficult to have much confidence in the form of the decline 
in the turnoff area around 2 Gyr in Figure The full 


isochrones at and at two rotation rates, the open stars 
in Figure show that the disappearance of an extended 
MSTO at ^2 Gyr may be even more abrupt than our 
Z = 0.008 extrapolation indicates. 

Finally, we suggest that the resolution of the disagree¬ 
ment between the range of ages suggested by the ex¬ 
tended MSTOs and by the observed spread in the sub¬ 
giant branch will require a detailed study. A range of 
rotation rates at constant age and composition will pro¬ 
duce a spread in masses at a fixed point in the turnoff 
region, just like a spread in ages. The mass ranges in 
these two scenarios are likely to be comparable and as 
large as ~10%. The extent to which rotation can reduce 
the spread in turnoff masses depends on the detailed be¬ 
havior of the stellar models with rotation and on the dis¬ 
tribution of initial rotation rates. Both deserve a much 
more detailed analysis than that presented here. 


The authors thank an anonymous referee for helpful 
suggestions, particularly on the expected distribution of 
projected rotational velocities. This work was performed 
in part under contract with the Jet Propulsion Labora¬ 
tory (JPL) funded by NASA through the Sagan Fellow¬ 
ship Program executed by the NASA Exoplanet Science 
Institute. 
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